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We give a summary of existing results on higher twist distribution amplitudes 
of vector mesons in QCD. Special attention is payed to meson mass corrections 
which turn out to be large. A “shopping list” is presented of most important 
nonperturbative parameters which enter distribution amplitudes. 


1 General framework 

The notion of distribution amplitudes refers to momentum fraction distribu¬ 
tions of partons in the meson in a particular Fock state with fixed number of 
components. For the minimal number of constituents, the distribution ampli¬ 
tude 4> is related to the Bethe-Salpeter wave function (jiBS by 


\k±\<fj. 



d^kj_ (l)Bsix,k±). 


( 1 ) 


The stajadard approach to distribution amplitudes, which is due to Brodsky and 
Lepageu, considers the hadron’s parton decomposition in the infinite momen¬ 
tum frame. A conceptually d^erent, but mathematically equivalent formalism 
is the light-cone quantizationa. Either way, power suppressed contributions to 
exclusive processes in QCD, which are commonly referred to as higher twist 
corrections, are thought to originate from three different sources: 

• contributions of “bad” components in the wave function and in particular 
of those with “wrong” spin projection; 

• contributions of transverse motion of quarks (antiquarks) in the leading 
twist components, given for instance by integrals as above with additional 
factors of 
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• contributions of higher Fock states, with additional gluons and/or quark- 
antiquark pairs. 

We take a somewhat different point of view and define light-cone distribution 
amplitudes as meson-to-vacuum transition matrix elements of nonlocal gauge- 
invariant light-cone operators. This formalism is convenient for the study 
of higher-twist distributions thanks to its gauge and Lorentz invariance and 
allows to solve all equations of motion explicitly, relating different higher-twist 
distributions to each other. We will find that all dynamical degrees of freedom 
are those describing contributions of higher Fock states, while all other higher- 
twist effects are given in terms of the latter without any free parameters. 

The report is divided into three sections, the first of which is introductory 
and the last two present the summary of distribution amplitudes up to twist 4. 
The expressions collected^ these sections are principally the result of recent 
studies reported in Refs, □’□tj which considerably extend the earlier analysis in 
□. We use a simplified version of the set of twist-4 distributions derived 
in B, taking into account only contributions of the lowest conformal partial- 
waves, and for consistency discard contributions of higher partial-waves in 
twist-3 distributions in cases where they enter physical amplitudes multiplied 
by additional powers of mp. Four-particle distributions of twist 4 start with 
higher conformal spin and must be put to zero to the present accuracy. The 
SU(3)-breaking effects are taken into account in leading-twist distributions and 
partially for twist-3, but neglected for twist-4. Explicit expressions are given 
for a (charged) p-meson. Distribution amplitudes for other vector-mesons are 
obtained by trivial substitutions. 

Throughout this report we denote the meson momentum by Pp and intro¬ 
duce the light-like vectors p and 2 ; such that 

1 rrii 

= —■ ( 2 ) 

The meson polarization vector is decomposed in projections onto the two 
light-like vectors and the orthogonal plane as 


ef,") = 


4A) . 


pz 


Pm 


2pz 


-'-Lm- 


( 3 ) 


We use the standard Bjorken-Drell convention!^ for the metric and the Dirac 
matrices; in particular 75 = and the Levi-Civita tensor is 

defined as the totally antisymmetric tensor with £0123 = 1. The covariant 
derivative is defined as Dp = l5p = dp — igAp, and we also use the notation 
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= d ^ + igAfj_ in later sections. The dual gluon field strength tensor is 
defined as 

1.1 Conformal partial wave expansion 

Conformal partial wave expansion in QCD i’i00’0lii parallels the partial 
wave expansion of wave functions in standard quantum mechanics, which al¬ 
lows to separate the dependence on angular coordinates from radial ones. The 
basic idea is to write down distribution amplitudes as a sum of contributions 
from different conformal spins. For a given spin, the dependence on the mo¬ 
mentum fractions is fixed by the symmetry. To specify the function, one has 
to fix the coefficients in this expansion at some scale; conformal invariance of 
the QCD Lagrangian then guarantees that there is no mixing between con¬ 
tributions of different spin to leading logarithmic accuracy. For leading twist 
distributions the mixing matrix becomes diagonal in the conformal basis and 
the anomalous dimensions are ordered with spin. Thus, only the first few “har¬ 
monics” contribute at sufficiently large scales (for sufficiently hard processes). 

For higher twist distributions, the use of the conformal basis offers a cru¬ 
cial advantage of “diagonalizing” the equations of motion: since conformal 
transformations commute with the QCD equations of motion, the correspond¬ 
ing constraints can be solved order by order in the conformal expansion. Note 
that relations between different distributions obtained in this way are exact: 
despite the fact that conformal symmetry is broken by quantum corrections, 
equations of motion are not renormalized and remain the same as in free theory. 

The general procedure to construct the, can formal expansion for arbitrary 
multiparticle distributions was developed inllJ’EJ. To this end each constituent 
field has to be decomposed (using projection operators, if necessary) in com¬ 
ponents with fixed (Lorentz) spin projection onto the light-cone. 

Each such component corresponds to a so-called quasiprimary field in the 
language of conformal held theories, and has conformal spin 

3 = \{l + s), (4) 

where I is the canonical dimension and s the (Lorentz) spin projection.^ 

Multi-particle states built of quasiprimary helds can be expanded in ir¬ 
reducible unitary representations with increasing conformal spin. An explicit 

= 3/2 for quarks and I = 2 for gluons; the quark field is decomposed as Tp = + pj— = 

(1/2) + (1/2)‘^^xp with spin projections s = +1/2 and s = —1/2, respectively. For the 

gluon field strength there are three possibilities: corresponds to s = +1, Pij.G^± has 

s = —1 and G_l_l, z^p^Gf^u correspond both to s = 0. 
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expression for the distribution amplitude of a multi-particle state with the low¬ 
est conformal spin j = ji + ■ • ■ + jm built of m primary fields with the spins 
jk is 



4^as (ai, q; 2, ..., a m 


Here ak are the corresponding momentum fractions. This state is nonde¬ 
generate and cannot mix with other states because of conformal symmetry. 
Multi-particle irreducible representations with higher spin j + n,n = 1,2,..., 
are given by polynomials of m variables (with the constraint Q^fe = 1 ) 

which are orthogonal over the weight function (i- 

1.2 Equations of motion 

We collect..here exact operator identities which can be derived using the ap¬ 
proach ofo and which present a nonlocal equivalent of the equations of motion 
for Wilson local operators. Taking suitable matrix elements one derives a set 
of relations between distribution amplitudes, which generally allow to express 
two-particle distributions of higher twist in terms of three-particle distribu¬ 
tions. The corresponding relations are taken into account in the results given 
in later sections. The rationale for presenting the operator relations themselves 
is that we found them useful in practical calculations: 



( 6 ) 


/. 


1 


d^u(x)ani^d(—x) = —z —— u(x)d(—x) + 

UXi/ 


-1 


dv vu{x)x^gGpv{vx)d{—x) 


— i dvu{x)x'^ 


L 


-1 


dv u{x)x'^gGpp{vx)a^^d{—x), 


( 8 ) 



/ 


1 


— I 


dv vu{x)x^gGpfj,{vx)a^j^d{—x), 


( 9 ) 


-1 


u{x)d{—x) — u{0)d{0) = dv u{tx)x^aa^x^gGp,i 3 {vx)d{—tx) 

t 
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( 10 ) 


+ i [ dtd°‘ {u{tx)aa0X^d{—tx)} . 

Jo 

In all cases gauge factors are implied in between the constituent fields, 

[x,y] =Pexp[ig f dt {x - y)f,A>^{tx + {1 - t)y)], (11) 

Jo 

and we introduced a shorthand notation for the derivative over the total trans¬ 
lation: 

da{u(tx)Td{-tx)} = -^{u(tx + y)Td{-tx + y)} (12) 

dy^ y^O 

with the generic Dirac matrix structure F. For simplicity, we omit operators 
involving quark masses, seeau. 

Two more relations are: 

d _ 

u{x)jij,d{—x) = dt —— u{tx)Jd{—tx) 

Jo OXfj, 

— I I ^^dGfj.u{vx)x^j'^d{—tx) 

+ J dtt J dvu{tx)^x'^gGfj,^{vx)j‘'j5 - x^x''gG^p{vx)Yl5'^d{-tx) 

-iepuais [ dttx''d°‘ [u{tx)^^^ 5 d{-tx)] , 

Jo 


(13) 


'0 

and similarly with an additional 75 : 


r 9 _ 

u{x)jnj5d{-x) = dt -— u{tx)J^od{-tx) 

Jo 

-I- J dtt J dvu{tx)^x'^gGij,,y{vx)'f'' — XpX''gG^p{vx)'f'^'^d{—tx) 

— i^tivag [ dttx''d°‘ [^{tx)^^d{—tx)] . 

Jo 


( 14 ) 


Finally, the following formula is sometimes useful: 


dXadx^^ 


u(a;)F d{—x) = —d^u{x)T d{—x) + u{x)[TaG + aGT]d{—x) 
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a /■! _ - 

— 2ix '^—— / dvvu{x)TGvfi{vx)d{—x)—2ix'^d^l dvu{x)TGi,p,{vx)d{—x) 
OX^J_l 7—1 

+ 2y J dt {1 + vt)u{x)Tx^x'^Gfj,p{vx)G^^(tx)d{—x) 

+ ix‘' J dv {1 + v^)u{x)r[Dp,G'^^^]{vx)d{—x) (15) 

where [Dp,G^^] = assuming summation over light flavors ijj. 


1.3 Meson mass corrections 

The structure of meson mass corrections in exclusive processes is in general 
more complicated than of target mass corrections in deep inelastic sc^tering 
in which case they can be resummed using the Nachtmann variable ta. For 
illustration, consider the simplest matrix element 


{0\u{x)^d{—x)\p {P, X)) = 
= fpmp(e^^'>x) [ 


(j){u) + —$(■«) + 0(a;"‘) 


(16) 


We assume that x^ but nonzero, is the twist-2 distribution 

amplitude and $(u) describes higher-twist corrections in which we want to 
calculate “kinematic” contributions due to nonzero p-meson mass. 

A common wisdom tells that meson mass corrections are related to con¬ 
tributions of leading twist operators. Indeed, conditions of symmetry and zero 
traces for twist-2 local operators imply 


(0| w^(j £) x)”d Ip (P, A)) = 
L J tw. 2 


= fpmpie^^^x) 


(17) 

4 n -I- 1 


where [.. .]tw .2 denotes taking the leading twist part (subtraction of traces, in 
this case) and ((On)) is the reduced matrix element related to the n-th moment 
of the leading twist distribution 

= f'du (2u - ir^iu) = ((On)). (18) 

Jo 

Expanding ( |l^ at short distances x ^ 0 and comparing with ( ITtI) , we find 
that the same reduced matrix element gives a contribution to the twist 4 
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distribution amplitude 


M(^) ^ fdu {2u - !)"$(«) = (19) 

Jo n + 3 ^ 

which is the direct analogue of Nachtmann’s correction for deep inelastic scat¬ 
tering. 

The result in jl^ ) is, however, incomplete. The reason is that in exclusive 
processes one has to take into account higher-twist operators containing full 
derivatives, and vacuum-to-meson matrix elements of such operators reduce, 
in certain cases, to powers of the meson mass times re 
of leading twist operators. For the case at hand, writf 


need matrix elements 


u(x)J:d{—x) = u{x)J:d{—x) 


u{x)Jd{—x) 


- tw.2 


X 

T 

T 


dt 




dxndx°‘ 


u{tx)Jd{—tx) -I- O(x^) 


dtr [u{tx)Jd[—tx)] 


contributions of operators with gluons -I- 0{x'^) 


( 20 ) 


where we used Eq. (^) to arrive at the last line. In the matrix element one can 
make the substitution —> —rn^. Expanding, again, at short distances, and 
comparing with the similar expansion of (0) we get an additional contribution 
to so that the corrected version of (H) becomes 


+ ((On))] + gluons. 


( 21 ) 


Assuming the asymtotic form of the leading-twist distribution amplitude (j), 
(j){u) = 6m( 1 — u), so that ((On)) = 3/[2(n -I- l)(n -|- 3)], this equation for 
moments is easily solved and gives 


$(w) =30 u^(1-m)^ 


2 2,4 2 /- 

5 Wp + ^^pC4 


( 22 ) 


where we have included the “genuine” twist 4 correction (term in ^ 4 ) due 
to the twist 4 quark-gluon oper^or, see definition in Eq. (^. The QCD 
sum rule estimate is ((4 ~ 0.15 113 , so that the meson mass effect onto the 
twist 4 distribution function is by a factor two larger than the “genuine” twist 4 
correction. This is an important difference to deep inelastic scattering, where 
the target mass corrections are small. 

The present discussion is still oversimplified and does not provide with 
a complete separation of meson mass effects. The major complication arises 
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( 23 ) 


because of contributions of operators of the type 


■ ■ ■ iiDf,Jd 


tw.2 


Such operators can be expressed in terms of operators with extra gluon fields, 
which means, conversely, that certain combinations of qGq operators reduce to 
divergences of leading twist operators and give rize to extra meson mass cor¬ 
rection terms. The corresponding corrections to twist 4 distributions involve, 
however, higher-order contributions in the conformal expansion of the distri¬ 
bution amplitudes of leading twist and do not affect the result in ([^), which is 
to leading conformal spin accuracy. A detailed analysis will be presented inU. 


2 Summary of chiral-even distributions 


Two-particle quark-antiquark distribution amplitudes are defined as matrix 
elements of non-local operators on the light-cone □: 


{0\u{z)-ffj,d(-z)\p {P,X))=fpmp 


Pf^- 




+e 


W 


P'Z Jq 

{V) / • 2 \ 1 ^ ■ z 

9x 


..2 


du {u, 


(24) 


and 


{{)\u{z)-ip'y^d{-z)\p (P,A)) = 


= o /p-/, 




■rud 


f duP^P'^gP{u,p'^).{25) 
Jo 


For brevity, here and below we do not show the gauge factors between the 
quark and the antiquark fields and use the short-hand notation 


^ = M — (1 — u) = 2u — 1. 

The vector and tensor decay constants fp and /J are defined, as usual, as 

(0|u(0)7^d(0)|p"(P, A)) = fpmpe^^\ (26) 

(0|h(0)a^,d(0)|p-(P, A)) = - e\^^Pp). (27) 

The distribution amplitude (/>|| is of twist-2, g^ and g^p are twist-3 and gs is 
twist-4. All four functions </> = {?!>||, 53 } are normalized as 









which can be checked by comparing the two sides of the defining equations 
in the limit > 0 and using the equations of motion. We keep the (tiny) 
corrections proportional to the u and d quark-masses niu and md to indicate 
the SU(3)-breaking corrections for K*- and <()-mesons. 

In addition, we have to define three-particle distributions: 


{Q\u{z)gG^v{vz)-ialbd{-z)\p (P, A)) = fpmpPa[pv 

. z 

pz 

■ z 

{0\u{z)gGp,,{vz)i-fad{-z)\p~ (P)) = fprUpPalp^e^^l 

+ fp^p—Z^\Pp9o.^ - 

pz 

+ fp^p Po.[PpZ^ - 


e-^ll-Ppe-^±l\Av,pz) 

Pv9ip\^{v,pz) 
-PuZp\^{v,pz), (29) 
-Ppe'll\^{'v,pz) 

Pvgip\Hv,pz) 

-pyZp\^{v,pz), (30) 


where 

A{v,pz) = J Pae-*P"(““-“‘^+"“»U(a), (31) 

etc., and a is the set of three momentum fractions a = {od, Ou, Qfg}. The 
integration measure is defined as 

J Va = J dad J dau J dag 5(1 — E“*)- (32) 

The distribution amplitudes V and A are of twist-3, while the rest is twist-4 
and we have not shown further Lorentz structures corresponding to twist-5 
contributions^. 

Calculation of exclusive amplitudes involving a large momentum-transfer 
reduces to evaluation of meson-to-vacuum transition matrix elements of non¬ 
local operators, which can be expanded in powers of the deviation from the 
light-cone. To twist-4 accuracy one can use the expression for the axial-vector 
matrix element in ( p5|) as it stands, replacing the light-cone vector Zp by the 
actual quark-antiquark separation Xp. For the vector operator, the light-cone 
expansion to twist-4 accuracy reads: 


{0\u{x)"fpd{-x)\p (P,A)} = 

*’We use a different normalization of three-particle twist-3 distributions compared toE 
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= fpmp 


Px 




TTl^X^ 

H— ^Hu,d) 




Px 


+ e\^ - P, 




1 e^^^x 

~ 2^^1p^ 


TOpy due’'^^^C{u, fi) 


where 


(33) 

(34) 


C(m) = gsiu) + (j)n{u) -2g^l\u), 

and A(it) can be related to integrals of three-particle distributions using the 
equations of motion. All distribution functions in (|^ ) are assumed to be 
normalized at the scale g? ^ x~^ (to leading-logarithmic accuracy). 

For the leading twist-2 distribution amplitude i^n we use 


^11 (u) = 6uu 


l + 3ak + 4-(5e-l) 


( 35 ) 


The parameters 2 defined as the local matrix elements 


)M^(i n z)dlp (P,A)) = (e(^)z)(pz)/pTOp-af, 


(0|u^(i D zydjp (P, A)) = z){pzYfpmp 17 + ^ 


1 12 


( 36 ) 


The numerical values are specified in Table |^. The expressions for higher-twist 
distributions given below correspond to the simpl^t self-consistent approxi¬ 
mation that satisfies the QCD equations of motionQ’Q: 


Three-particle distributions of twist-3: 


V(a) = (ttd - au)adauag, 

A(a) = 360C3ada«ag[l+w^^(7ag-3)]. 


(37) 

(38) 


Two-particle distributions of twist-3: 


g^^(u) = 6 mm 


l-l-a)'^-|-<i Tal+oCsl “ 1) 


16 


16 


- 1-6 5 + ( 3 mm -|- mI n m -I- mI n m) -I- 6 5 _ (Mlnu — MlnM), 


(39) 
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9S:’(«) = j(i+?“)+«! +(f 4+ 5 C 3 ) (se - 1 ) 


112 


15 


C3( 


(3-30^2 + 35^^) 


3 ~ 3 ~ 

+ - (2 -h In u -h In u) -h - S- (2^ + Inu — Inu). 


(40) 


• Three-particle distributions of twist-4: 
$(a) = 

$(a) = 






30 (l-ag)ag, 


30(au - ad)al 


[2 1 1 

^(a) = [ 3 C 3 + gC4j 120 a„adaff, 


4-(a) = 0 . 

• Two-particle distributions of twist-4: 
A(u) = 


(41) 


4 , 4 II 20 8 

5 + 1^ “2 + yC4 + gC3 


gsiu) = 6 u{l - u) + 
C(u) = 


1 II 10 ^ 20 ^ 

y4 + yC4-yC3 


3 2 II 10 ^ 10 ^ 

2-7«I + yC4 + yC3 


30m 2(1 -^t)^ 

(i-sa, 

(1-3^^), 


(42) 


where the dimensionless couplings ^3 and (^4 are defined as local matrix ele¬ 
ments 


{0\ugGfj,„'yal5d\p (P, A)) = fpiripCs 


A>') (p p - \ 

a J 14 2 nf^p yai4 j 


{PaPp. - 




eWn -e^^'>a 

yotu ^ 1 / yafi 


,( 43 ) 


and have been estimated from QCD sum- 


rules 00 . ojY is defined as 


'\u^(gGa 0 Z°‘i(D z) - (i D z)gGai 5 Z°‘)d\p (P, A)) = 


= iipz) eYpmpfp— CaujY + O{z 0 ) (44) 
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Figure 1: The twist-2 distribution amplitude 0|| (u) for p, K* and (p mesons. Renormalization 
point is p = 1 GeV. 



and is defined as 


(0|u^75 


iDz,gGfj,^z^ 


d\p {P, A)) = -{pzYe\lmpfpC,3 ( - + —Wg ) +0{zf}) 


(45) 

The terms in (5± and i5± specify quark-mass corrections in twist-3 distributions 
induced by the equations of motion. The numerical values of these and other 
coefficients are listed in Tables and s Note that we neglect SU(3)-breaking 
effects in twist-4 distributions and in gluonic parts of twist-3 distributions. In 
Fig. 1 we plot the leading-twist distribution amplitude for p, K* and (j) 
mesons. 


3 Summary of chiral-odd distributions 

There exist four different two-particle chiral-odd distributions S defined as 
{0\u{z)ap„d{-z)\p-{P,X)) =ifj - e^^Pp) [ duP^P'^(t)±{u, pY 




+ f duP^P'^h3{u,pY 

2 P P- z Jo 


,(46) 


*^10 the notation of Ref.a, = o;^, = ^3, and = (3/2S). 
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Table 1: Masses and couplings of vector-meson distribution amplitudes, including SU(3)- 
breaking. In cases where two values are given, the upper one corresponds to the scale 
= 1 GeV^ and the lower one to = 5 GeV^, respectively. We use ms{l GeV) = 150 MeV 
and put the u and d quark mass to zero. 


V 



^u,d 

<t> 

/y[MeV] 

198 ±7 

226 ± 28 

226 ± 28 

254 ±3 

/^[MeV] 

160 ± 10 

185 ± 10 

185 ± 10 

215 ± 15 

152 ±9 

175 ±9 

175 ±9 

204 ± 14 

II 


0.19 ±0.05 

-0.19 ±0.05 


af 

0 

0.17 ±0.04 

-0.17 ±0.04 

0 

II 

0.18 ±0.10 

0.06 ±0.06 

0.06 ±0.06 



0.16 ±0.09 

0.05 ±0.05 

0.05 ±0.05 

0 ±0.1 

aj- 


0.20 ±0.05 

-0.20 ±0.05 


0 

0.18 ±0.05 

-0.18 ±0.05 

0 


0-2 

0.20 ±0.10 

0.04 ±0.04 

0.04 ±0.04 

0 ±0.1 


0.17 ±0.09 

0.03 ±0.03 

0.03 ±0.03 




0.24 

0.24 

0.46 

s+ 

0 

0.22 

0.22 

0.41 



-0.24 

0.24 


5- 

0 

-0.22 

0.22 

0 



0.16 

0.16 

0.33 

s+ 

0 

0.13 

0.13 

0.27 



-0.16 

0.16 


6 - 

0 

-0.13 

0.13 

0 
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Table 2: Couplings for twist-3 and 4 distribution amplitudes for which we do not include 
SU(3)-breaking. Renormalization scale as in the previous table. 



Cs 

/ A 

0)3 

0)3^ 

T 

^3 

C4 

a 

CJ 

V 

0.032 

-2.1 

3.8 

7.0 

0.15 

0.10 

-0.10 


0.023 

-1.8 

3.7 

7.5 

0.13 

0.07 

-0.07 


{Q\u{z)d{-z)\p {P,\)) = 

(e(A) (47) 

\ '^P J Jo 

The distribution amplitude (j)± is twist-2, and are twist-3 and /13 is 
twist-4. All four functions (j) = {4>±, ^ 3 } sre normalized to 

/ du (j){u) = 1. 

Jo 

Three-particle chiral-odd distributions are defined to twist-4 accuracy as 
{0\u{z)aapgGf,„{vz)d{-z)\p~{P, A)) = 

g(A) . 2 

= /J 2(p ■ z) ~ PdPl^sL - PaPt^dp^ + (v, pz) 

+ ppe^llg^^]T^‘^\v,pz) 

+ fjml[pf,e^llgj^ - p^e^l^g^^ - p^e^llgj^ + p^e^llg^^]TP {v,pz) 
fp^l r„ „(A), (A), ^ ^ (A)^ , ^(A), 1„(4). 


I •'P ^P\ (A) (A) {A) Trrii^l)/■ \ 

+ [PaPpe)_^z^ - ppPf,e\^z^ - PaPue^pZf, + PpPue\^z^\T^ {v,pz) 


fT'rn^ 

+ — — -[PaPpe^llzp -ppp^e^^lza - PaPve^^lzp + ppPi,e^^lza]Tl^\v,pz) 

( 48 ) 


AA). 


AA). 


JA). 

AA). 


AA) T.rp(4:)i 


and 

{Q\u{z)gG^^{vz)d{-z)\p~{P,X)) = ifjmlle^^^p^ - e^^lpfj,]S{v,pz), 
{0\u{z)igGf,^{vz)-f5d{-z)\p-{P, \)} = - e^llp^]S{v,pz).{A9) 

Of these seven amplitudes, T is twist-3 and the other six are twist-4. 
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The light-cone expansion of the non-local tensor operator can be written 
to twist-4 accuracy as 

{Q\u{x)a^ud{-x)\p~{P,X)) = 



where Bt and Ct are expressed in terms of the distribution amplitudes defined 
above as 

Bt(m) = h\^\u) - - ^hsiu), 

Ct{u) = hsi^u) - (51) 


and At can be related to integrals over three-particle distribution functions 
using the equations of motion. 

For the leading twist-2 distribution amplitude (p± we use 


(j)±{u) = 6 uu 


1 + Sa]*" C + 



(52) 


with parameter values as specified in Table Q the definitions of aj ~2 are analagu- 
ous to Eqs. (^). The expressions for higher-twist distributions given below 
correspond to the simDlest self-consistent approximation that satisfies all QCD 
equations of motion Qa: 


• Three-particle distribution of twist-3: 

T(a) = 540 ( 3 UJ 3 (ad - au)adaual. (53) 


• Two-particle distributions of twist-3: 


hi^\u) = 6 uu 


1 -l- ctj*"^ -)- 


1 


■C3W3 


(5C^ - 1) 


-I- 315+ (3uu -I- uInrt -I- uInu) -I- 315_ (itInu — ulnw), (54) 
(u) = 3^2 -b ^ ^(3^2 - 1) + (5^^ - 3) 
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(55) 


1 ^ Q 

+ ~ + 2 '^+ ^ Inw/u) 

3 

+ - (5_ ^ (2 + Inw + Inw). 

• Three-particle distributions of twist-4: 

= Tt\a) = 0 , 

T^^\a) = SOCliaa - au)al 
Tt\a) = -30(1 {ad - ctu)al, 

S{a) = 30(Jil-ag)al, 

S{a) = 30CJ(1 - Q!g)Q!g. 


• Two-particle distributions of twist-4: 


hsiu) 


At{u) 


6u(l — u) -I- 5 

30it2(l - uf 


4 W S4 yg “2 

2 4 _L 4^-r 

5 + + 3^4 


(1 - 3e^), 



(56) 


(57) 


The constant wj is defined as 

{0\ua^^z‘'{gG^^Zf3{i D z) - {i D z)gG^^Zjd)d\p~{P,\)) = 

= {pzf{e^^h)mlf (58) 

The constants (J' and (J' are defined as 

(0|%G^.d|p-(P, A)} = - ei^^P^), 

{0\ugGfj,Jj5d\p~{P,X)) = ifpinlCJiel^^P^ - ei^'>Pp). (59) 

Numerical values for all parameters are given in Table As in the chiral-even 
case, we neglect SU(3)-breaking corrections in twist-4 distributions. In Fig. 2, 
we plot the leading-twist distribution amplitude (j)± for p, K* and (j) mesons. 
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Figure 2: Twist-2 distribution amplitude = 1 GeV) for p, K* and <j) mesons. 
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